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Abstract 
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On 

ON \ We report experiments on thermally driven convection in an inclined layer of large 

aspect ratio in a fluid of Prandtl number a « 1. We observed a number of new 
qj • nonlinear, mostly spatio-temporally chaotic, states. At small angles of inclination 

Q . we found longitudinal rolls, subharmonic oscillations, Busse oscillations, undulation 

chaos, and crawling rolls. At larger angles, in the vicinity of the transition from 
CN | buoyancy- to shear-driven instability, we observed drifting transverse rolls, localized 

bursts, and drifting bimodals. For angles past vertical, when heated from above, we 
found drifting transverse rolls and switching diamond panes. 
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Rayleigh-Benard convection (RBC) of a vertical fluid layer heated from below has long served 
as a paradigm for pattern forming systems Variations that alter the symmetries, such as ro- 
tation around a vertical axis [TJQ] and vertical vibrations [|5p continue to lead to important insights. 
Another case, of particular meteorological and oceanographic interest, is RBC of a fluid layer in- 
clined with respect to gravity. This system is not only well suited for the study of buoyancy and 
shear flow driven instabilities, but may also serve, along with liquid crystal convection [@], as a 
paradigm for anisotropic pattern forming systems. 

As with RBC, the onset of inclined layer convection (ILC) occurs when the temperature dif- 
ference, AT, across the layer is sufficient for convection rolls to form. The main difference from 
RBC is that, in ILC, the patternless base state is characterized not only by a linear temperature 
gradient but also by a symmetry-breaking shear flow. As shown in Fig. [j], the component of grav- 
ity tangential to the fluid layer, gy, causes buoyant fluid to flow up along the warm plate and down 
along the cold plate. For small angles of inclination 6, buoyancy dominates over shear flow (g^ 
large, gn small), and the primary instability is to longitudinal rolls (LR) whose axes are aligned 
with the shear flow direction [[7p. With increasing 9, buoyancy effects decrease and for 9 > 90° 
buoyancy is stabilizing. Above a critical angle 9 C (< 90°) the shear flow causes a primary insta- 
bility to transverse rolls (TR) with roll-axes perpendicular to the shear flow [|7p. The few prior 
experiments (TDQ on ILC showed reasonable agreement with the linear theory p7|JTT]— [T3t] . These 
experiments also demonstrated that LR are unstable to some form of undulations [R-ITT1], in quali- 



tative agreement with theory [ |l4|13| ], but the quantitative details of the state were inaccessible due 
to experimental limitations. 

Here we report the first experimental results on pattern formation in ILC for large aspect- 
ratio systems in a range of inclination angles 0° < 9 < 120°, i.e., from horizontal (heated from 
below) to past vertical (heated from above). We found many unpredicted states when increasing 
AT above the critical temperature difference. For 0° < 9 < 77.5° we observed longitudinal 
rolls, subharmonic oscillations, Busse oscillations, undulation chaos, and crawling rolls. In the 
neighborhood of the codimension two point [ fl~T| ] for thermal and shear driven instability (77.5° < 
9 < 84°), we observed drifting bimodals, drifting transverse rolls, and localized longitudinal and 
transverse bursts. For inclinations 9 > 84° we found drifting transverse rolls, switching diamond- 
panes, and longitudinal bursts. Most of these novel states were spatio-temporally chaotic and were 
found very close to onset, where theoretical progress should be possible. 

Experiment: Our experimental apparatus consisted of a water-cooled pressure chamber con- 
taining a convection cell of diameter 10 cm, subdivided into two large aspect ratio rectangular 
cells. The experimental design was similar to the one described in [[Rp. The optically flat upper 
and lower plate of the convection cell consisted of 1 cm thick single crystal sapphire and single 
crystal silicon, respectively. The sapphire plate was cooled by a water bath, while the silicon plate 
was heated by an electric film heater. The convection patterns were visualized by the usual shad- 
owgraph technique JT4t ] . The sidewalls were constructed of 9 layers of notebook paper, providing 
the best possible thermal matching between cell boundaries and the fluid. As measured interfer- 
ometrically, the plates were parallel to ±0.5 fxm. The convection cell was housed in a pressure 
chamber, which held both the cooling water and the convecting gas to (41.37± 0.01) bar, regulated 
to ±5 x 10~ 3 bar. The temperatures of the two plates were regulated to ±0.0003°C. Throughout 
the experiment the mean temperature was kept constant at (27.00 ± 0.05)°C. We determined the 
cell height d by measuring the pattern wavenumber at onset for 9 < 60° and comparing it with the 
theoretical value of q c = 3.117<i. We found d = (710 ± 2)fim and d = (702 ± 2)//m for two sets 
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of experiments. The two convection cells had a size Ti ~ (21 x A2)d 2 and T 2 — (14 x A8)d 2 . For 
all data, the Prandtl number was a = v / ' k = 1.07 as determined from [|T4|], with the kinematic 
viscosity v and thermal diffusivity k. The vertical thermal diffusion time was r v = d 2 / k = 3.0 s. 
Inclines from 0° (horizontal) to 120° (30° past vertical) were possible, with an accuracy of ±0.02°. 
Following [Q] we calculated the Boussinesq number V(9) for the corresponding horizontal layer 
to estimate non-Boussinesq effects. At AT C (#) for 9 < 70° we found V{9) < 1.0, putting the 
flow into the Boussinesq regime. For larger angles V increased linearly to 3.0 for the largest tem- 
perature differences investigated. We observed the same convective patterns in both convection 
cells. 

Onset of convection: In ILC, the forward bifurcation to LR is predicted to occur at the critical 
Rayleigh number R c (9) = R t c (0°)/cos9 where i?*(0°) = 1708 = ad 3 gAT c /niy (a is the thermal 
expansion coefficient and AT C is the critical temperature difference). The threshold for the forward 
bifurcation to shear driven TR at large inclination angle is more complicated, and can only be 
determined numerically [ |I^|Dl , |i5| ]. 

We determined AT C for convection by quasi- statically increasing AT in steps of 1 mK every 
20 minutes past the point where convection was observable and then decreasing the temperature 
difference similarly. For all angles we observed forward bifurcations. Figure || shows the measured 
R c (9), as well as the theoretically predicted onsets for both the buoyancy-driven (longitudinal) and 
the shear-driven (transverse) instabilities [[13]]. We found agreement with theory for the onset of 
LR: the experimentally observed value was R c (9) = R e c (0°)/cos9 with R°(0°) = 1687 ± 24. 
We did not, however, observe the theoretically predicted stationary TR, but instead drifting TR 
(DTR) at a slightly larger critical Rayleigh number. The drift down the incline may be attributed 
to the broken symmetry across the layer which is caused by the temperature dependence of the 
fluid parameters (non-Boussinesq effects). Very interesting is the vicinity of the theoretically 
predicted codimension two point at 9 C = 77.76° JT3|] , where LR and TR have the same onset value. 
Experimentally, we found a forward bifurcation to DTR above 9 C = (77.5 ± 0.05)°, and in the 
range 77.5° < 9 < 84° DTR lost stability to drifting bimodals (DB) above e « 0.001. As shown 
in Fig. |, DB consist of a superposition of LR and DTR. Here e = (AT - AT c (0))/AT c (0) is the 
reduced control parameter. Theoretically, Fujimura and Kelly JTTt ] predicted a forward bifurcation 
to transverse rolls, which lose stability to bimodals at e « 0.001 in a narrow angular region. 
We find good agreement with these predictions, but with the difference that the experimentally 
observed patterns are drifting. 

Nonlinear states: Figure f| shows the measured phase boundaries for the ten observed nonlin- 
ear convective states. At low angles (9 < 13°), LR are stable up to e ~ 1 above which the novel 
state of subharmonic oscillations (SO) sets in. These oscillations are characterized by a pearl- 
necklace-like pattern of bright (cold) spots that travel along a standing wave pattern of wavy rolls. 
As shown in Fig. |5], these oscillations appear in patches whose location changes in time. Typi- 
cal frequencies of the oscillations were measured to be 1 to 3 cycles per r v . A recent theoretical 
analysis has shown agreement with this value [|T8|]. With further increase in e, localized patches 
of traveling oscillations burst intermittently. Within 0(r v ) the amplitude of the rolls' waviness 
increases, the pattern tears transverse to the rolls as shown in the upper left corner of Fig. BL and 
fades away leaving an almost parallel roll state. 

For 9 ~ 10° and e > 4, we observed patches of the well-known Busse oscillations (BO) 
coexisting with patches of the SO. As shown by the dotted line in Fig. |], our data for the onset 
of the BO agrees well with the theoretical prediction calculated for a = 0.7 JT2]]. It is surprising, 
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however, that both oscillations (SO and BO) coexist as localized patches in the same cell. 

At intermediate angles (25° < 9 < 70°), where the initial instability is to LR we found with 
increasing e that LR were unstable to undulations. Although the experimentally determined value 
for the instability e ~ 0.01 agrees well with the theoretical prediction (see Fig. 0) JL^JTIIJL^ ], 
we did not observe a stationary pattern of undulations, but a defect-turbulent state of undulation 
chaos (UC). A snapshot of UC is shown in Fig. ||a. At e > 0.11, the UC begins to "twitch" 
in the direction transverse to the rolls on time scales 0(r v ). With increasing e, the amplitude 
of the twitching increases and the rolls eventually tear, with the ends "crawling" in the direction 
transverse to the original rolls. A snapshot of crawling rolls (CR) is shown in Fig. |^b. 

In the vicinity of the codimension-two point, at 9 C , we observed drifting bimodals quite close 
to onset. As shown in Fig. f|, for small angles the existence region of the pure DB is limited by 
localized transverse bursts (TB), while for large angles by DTR. A snapshot of transverse bursts 
and the evolution of a single burst is shown in Fig. [7]. In this region of phase space the LR occur in 
patches that grow and decay intermittently while TB nucleate in high amplitude LR-regions. As 
shown in the time series in Fig. [7], TB grow over the period of a few r v and then decay rapidly. 
Above e w 0.8 the DB are unstable to localized longitudinal bursts (LB) as shown in Fig. [8|a. As 
shown in Fig. |8|b-j, a few longitudinal rolls grow locally to large amplitude and then quickly fade. 
With both types of bursts, the bursts increase both in density and frequency when e is increased, 
eventually developing into a turbulent state at e > 1. 

Past 90°, we continued to observe shear driven convection patterns. DTR are the primary 
instability; however, they are unstable to switching diamond panes (SDP) at e ~ 0.07. The state is 
characterized by spatio-temporally chaotic switching on time-scales of 0(r v ) from +45° to —45° 
of large amplitude regions of DTR, as seen in Fig. |9|a. At e > 0.1 SDP are unstable to LB, which 
in this region of phase space are denser but travel less distance than in TR, as shown in Fig. |^b. 

Conclusion: Inclined layer convection in the weakly nonlinear regime displays a rich phase 
diagram, with ten different states accessible over the range of parameters investigated. The phase 
space naturally divides into several regions of characteristic behavior which have so far been char- 
acterized semi-quantitatively. All states but longitudinal and transverse rolls are spatio-temporally 
chaotic. Most instabilities occurred very close to onset and further theoretical description should 
be possible. Especially interesting is the bursting behavior, which may be related to turbulent 
bursts in other shear flows []19[]. 
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theoretical descriptions of various states. E.B. acknowledges the kind hospitality of Prof H. Levine 
at the University of California at San Diego where part of this manuscript was prepared. We 
gratefully acknowledge support from the NSF under grant DMR-9705410. 
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FIG. 1. Schematic of the base flow, (a) Heated from below and (b) heated from above; cell height d, 
gravitational acceleration g, shear flow v, and temperature difference AT = T2 — T\ with T% > T\. 
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FIG. 2. Onset of longitudinal rolls (o) and drifting transverse rolls (A). Also plotted are the predicted 



onsets for longitudinal (dashed) and transverse rolls (solid) [15]. 
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FIG. 3. Digitally enhanced shadowgraph image of bimodals drifting from left to right in Cell 1 for 
6 = 77.6°, e = 0.01. The left side of the cell is higher than the right, with warm up-flow to the left and 
cold down-flow to the right. The rolls at the edges of the cell are caused by sidewall imperfections. MPEG 
movies of all states are available online [16]. 
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FIG. 4. (e, 6) phase-space showing the boundaries between the different nonlinear states. LR (longitu- 
dinal rolls), BO (Busse oscillations), SO (subharmonic oscillations), UC (undulation chaos), CR (crawling 
rolls), DTR (drifting transverse rolls), DB (drifting bimodals), LB (longitudinal bursts), TB (transverse 
bursts), and SDP (switching diamond panes). The dotted line is the predicted onset of Busse oscillations 
for a 
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the dashed line is the predicted onset of undulations [15], and the solid lines are guides 



to the eye. Open circles (UC) were measured via defect density Q17|], open diamonds (SDP) were measured 
via correlation length JT7[], and the remainder were measured visually. The inset shows a magnification of 
the codimension-two region. 
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FIG. 5. Contrast enhanced shadowgraph image of subharmonic oscillations in Cell 1, with a turbulent 
burst in the upper left corner. 6 = 17°, e = 1.5. Warm up-flow to the left and cold down-flow to the right. 
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FIG. 6. Digitally enhanced shadowgraph images of convection states at 9 = 40° in Cell 1 . (a) Undula- 
tion chaos at e = 0.07. (b) Crawling rolls at e = 0.88. Warm up-flow to the left and cold down-flow to the 
right. 
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FIG. 7. Digitally enhanced images of (a) transverse bursts in spatio-temporally chaotic longitudinal 
rolls, at 6 = 77° and e = 0.04 in Cell 1. (b - u) time-evolution of a single burst at time-intervals 0.18r v . 
Warm up-flow to the left and cold down-flow to the right. 
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FIG. 8. Digitally enhanced images of (a) longitudinal bursts at 6 = 79° and e = 0.10 in Cell 1. (b - 
j) time-evolution of a single burst at time-intervals 0.09r v . Warm up-flow to the left and cold down-flow to 
the right. 
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FIG. 9. Digitally enhanced shadowgraph images of (a) switching diamond panes (e = 0.1, 9 = 100°) 
and (b) longitudinal bursts within diamond panes (e = 0.19, 9 = 100°) in Cell 1. Warm up-flow to the left 
and cold down-flow to the right. 
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